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ABSTRACT 

We investigate a hierarchical structure formation scenario in which galaxy stellar cores are created 
from the binding energy liberated by shrinking supermassive black hole (SMBH) binaries. The binary 
orbital decay heats the surrounding stars, eroding a preexisting stellar cusp oc r^^. We follow the 
merger history of dark matter halos and associated SMBHs via cosmological Monte Carlo realizations 
of the merger hierarchy from early times to the present in a ACDM cosmology. Massive black holes 
get incorporated through a series of mergers into larger and larger halos, sink to the center owing to 
dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, and 
form a binary system. Stellar dynamical processes drive the binary to harden and eventually coalesce. 
A simple scheme is applied in which the loss cone is constantly refilled and a constant density core forms 
due to the ejection of stellar mass. We find that a model in which the effect of the hierarchy of SMBH 
interactions is cumulative and cores are preserved during galaxy mergers produces at the present epoch a 
correlation between the 'mass deficit' (the mass needed to bring a flat inner density profile to a cusp) 
and the mass of the nuclear SMBH, with a normalization and slope comparable to the observed relation. 
Models in which the mass displaced by the SMBH binary is replenished after every major galaxy merger 
appear instead to underestimate the mass deficit observed in 'core' galaxies. 

Subject headings: cosmology: theory - black hole physics - galaxies: evolution - quasars: general 



1. INTRODUCTION 

The strong link observed between the masses of super- 
massive black holes (SMBHs) residing at the center of 
most nearby galaxies and the gravitational potential wells 
that host them (Ferrarese & Merritt 2000; Gebhardt et al. 
2000; Ferrarese 2002) suggests a fundamental mechanism 
for assembling black holes (BHs) and forming spheroids 
in galaxy halos. In popular cold dark matter (CDM) 
'bottom-up' cosmogonies, the host galaxies of SMBHs ex- 
perience multiple mergers during their lifetime (Kauff- 
mann & Haehnelt 2000; Menou, Haiman, & Narayanan 
2001), with those between comparable-mass systems ('ma- 
jor mergers') expected to result in the formation of ellip- 
tical galaxies (Barnes 1988). If the hierarchical build-up 
of SMBHs traces far up in the dark halo merger tree, bi- 
nary BH systems may be widespread (Volonteri, Haardt, 
& Madau 2003, hereafter Paper I). As galaxies hosting 
nuclear BHs merge, the holes will sink to the center of 
the merger remnant owing to dynamical friction from field 
particles, and form a bound pair. 

If was first proposed by Ebisuzaki, Makino, & Okumura 
(1991) that the heating of the surrounding stars by a de- 
caying SMBH pair would create a low-density core out of a 
preexisting cuspy stellar profile. In a stellar background a 
'hard' binary shrinks by capturing the stars that pass close 
to the holes and ejecting them at much higher velocities, a 
super-elastic scattering process that depletes the nuclear 
region. Rapid coalescence eventually ensues due to the 
emission of gravitational radiation. Observationally, there 
is clear evidence in early-type galaxies for a systematically 
different distribution of surface brightness profiles, with 
faint ellipticals showing steep power-law profiles (cusps). 



while bright ellipticals have much shallower stellar cores 
(Lauer et al. 1995; Gebhardt et al. 1996; Faber et al. 1997; 
Ravindranath et al. 2001). 'Core' galaxies have high veloc- 
ity dispersions and exhibit a definite break in the bright- 
ness profile at some radius that scales with luminosity as 
oc Ly^^ (Faber et al. 1997). Detailed N-body simulations 
have confirmed the cusp-disruption effect of a hardening 
BH binary (Makino & Ebisuzaki 1996; Quinlan & Hern- 
quist 1997; Milosavljevic & Merritt 2001), but have shed 
little light on why bright ellipticals have lower central con- 
centrations than do faint ellipticals. 

The role of binaries in shaping the central structure of 
galaxies can be best understood within the framework of 
a detailed model for the hierarchical assembly of SMBHs 
over cosmic history (Paper I; Haehnelt & Kauffmann 
2002), particularly if the damage done to a stellar cusps 
by decaying BH pairs is cumulative and nuclear cores are 
preserved during galaxy mergers. In this paper we study 
the effects of hierarchical mergers of halo-|-SMBH systems 
on the inner density profiles of galaxies using the machin- 
ery for following the growth and dynamics of SMBHs de- 
veloped in Paper I. We show that stellar cusps can be 
efficiently destroyed over cosmic time by decaying SMBH 
binaries if stellar dynamical processes are able to shrink 
the binary down to a separation < 10% of the separation 
at which the binary becomes 'hard'. More massive ha- 
los have more massive nuclear BHs and experience more 
merging events than less massive galaxies: hence they suf- 
fer more from the eroding action of binary SMBHs and 
have larger cores. 

2. ASSEMBLY AND GROWTH OF SMBHs 



^Department of Astronomy & Astrophysics, University of California, Santa Cruz, CA 95064. 
■^Dipartimento di Scienze, Universita dellTnsubria/Sede di Como, Italy. 



2 



Formation of Galaxy Stellar Cores 



We briefly summarize here the main features of our sce- 
nario for the hierarchical growth of SMBHs in a ACDM 
cosmology (see Paper I for a thorough discussion). The 
merger history of 220 parent halos with present-day masses 
in the range 10^^ < Mq < 10^^ Mq is tracked back- 
wards with a Monte Carlo algorithm based on the ex- 
tended Press-Schechter formalism. Compared to Paper 
I, we use an improved version where the most massive ha- 
los are broken up into as many as 280,000 progenitors by 
z = 20. We adopt a two-component model for galaxy ha- 
los. The dark matter is distributed according to a NFW 
profile (Navarro, Frenk, & White 1997), 



PDM(r) 



M 



47rr(r + rvir/c)2/(c) 



(1) 



where r^iT is the virial radius, c is the halo concentration 
parameter and 



/(c)=ln(l + c) 



1 + c 



(2) 



Following Bullock et al. (2001), the mean concentration is 
assumed to scale with halo mass M and redshift of collapse 
z as 

__ 9 / M \-°-^" 

l + z\8xW^^hMQ) ■ ^' 

Individual halos have a log-normal distribution with dis- 
persion about the mean A (log c) = 0.18. During the 
merger of two halo+BH systems of comparable masses, 
dynamical friction against the dark matter background 
drags in the satellite hole towards the center of the newly 
merged system, leading to the formation of a bound BH 
binary in the violently relaxed stellar core. The dynamical 
friction timescale depends on the orbital parameters of the 
infalling satellite, which we take from van den Bosch et al. 
(1999). At late epochs, most of the BH pairs have unequal 
masses, with mass ratios ranging between 10% and 20% 
(Paper I). 

The subsequent evolution of the binary is determined 
by the initial central stellar distribution. We model this 
as a singular isothermal sphere (SIS) with one-dimensional 
velocity dispersion cr* and density 



P*{r) = 



27rGr2 ' 



(4) 



each major merger the more massive hole accretes at the 
Eddington rate a gas mass that scales with the fifth power 
of the circular velocity of the host halo, 



3.6 X 10'^ M0 IC 



(6) 



where V^ciso is the circular velocity of the merged system 
in units of 150 kms~^.^ The normalization is fixed a pos- 
teriori in order to reproduce the observed local tubh — 
relation (Ferrarese 2002). The present-day mass density of 
nuclear SMBHs accumulates mainly via gas accretion, with 
BH-BH mergers playing only a secondary role. This model 
was shown in Paper I to reproduce remarkably well the ob- 
served luminosity function of optically-selected quasars in 
the redshift range 1 < z < 5. 

If the merging timescales of SMBH binaries are longer 
than the characteristic timescale between major galaxy 
mergers, then interactions between the binary and new 
infalling holes will be likely. In Paper I we discussed a 
scenario where binaries decay efficiently both as a result 
of mass ejection from a cuspy stellar density profiles and, 
at very high redshifts, due to triple BH interactions. In 
the next section we expand upon Paper I and describe a 
scheme for generating low-density cores out of a preex- 
isting stellar cusps from the binding energy liberated by 
shrinking SMBH binaries. 

3. DYNAMICAL EVOLUTION OF SMBH BINARIES 

Consider a binary with BH masses m\ > m2 and semi- 
major axis a{t) in an isotropic background of stars of 
mass ^ TO2 and density p*{r). The binary will ini- 
tially shrink by dynamical friction from distant stars act- 
ing on each BH individually. But as the binary separation 
decays, the effectiveness of dynamical friction slowly de- 
clines because distant encounters perturb only the binary 
center of mass but not its semimajor axis. The BH pair 
then hardens via three-body interactions, i.e., by captur- 
ing and ejecting at much higher velocities the stars passing 
by within a distance ~ a ('gravitational slingshot'). The 
system becomes hard when the value of a falls below 



a-h = 



Gm2 
4^ 



(7) 



The stellar velocity dispersion is related to the halo cir- 
cular velocity 14 at the virial radius following Ferrarese 
(2002), 

logyc = (0.88 ± 0. 17) log cr*-F (0.47 ±0.35). (5) 

We truncate the stellar SIS at 0.16 rvir in order for the total 
stellar mass fraction to equal the universal baryon fraction, 
Qii/Qm- In our model pregalactic 'seed' holes form with 
intermediate masses (m, = 150 Mq) in (mini)halos col- 
lapsing at z = 20 from rare 3.5-tT peaks of the primordial 
density field (Paper I; Madau & Rees 2001). The assumed 
'bias' assures that almost all halos above 10^^ Mq actu- 
ally host a BH at all epochs. We found little change in the 
z < 5 results in a test model case with m, = 1000 Mq. In 

^In eq. (6) a typographical error that appeared in the coefficient of eq. 



(Quinlan 1996). We assume that the 'bottleneck' stages 
of the (bound) binary shrinking occur for separations 
a < tth', during a galactic merger, after a dynamical friction 
timescale, we place the BH pair at o/, and let it evolve. 

The hardening of the binary modifies the stellar den- 
sity profile, removing mass interior to the binary orbit, 
depleting the galaxy core of stars, and slowing down fur- 
ther hardening. If A^ej is the stellar mass ejected by the 
BH pair, the binary evolution and its effect on the galaxy 
core are determined by two dimensionless quantities: the 
hardening rate 

H=^iK (8) 



Gp* dt a 



and the mass ejection rate 



J = 



1 



dMe 



(mi -H TO2) rfln(l/a)' 
(13) in Paper I has been corrected. 
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A third dimensionless quantity, the eccentricity growth 
rate, was shown by Quinlan (1996) to be unimportant 
and will be neglected here. The quantities H and J 
can be found from scattering experiments that treat the 
star-binary encounters one at a time. Following Quinlan 
(1996), we take H = 15 independent of a (for a < ah) 
and of the binary mass ratio. The quantity J is instead 
a function of binary separation and has some dependence 
on mi/ 1712- We inter /extra-polate Quinlan's numerical re- 
sults to obtain J{a) for different mass ratios (cf. Paper 
I). The merger timescale is computed adopting a simple 
semi-analytical scheme that qualitatively reproduces the 
evolution observed in N-body simulations (Merritt 2000; 
Paper I). We assume that the stellar mass removal creates 
a core of radius Vc and constant density pc = p^{rc), so 
that the total mass ejected as the binary shrinks from 
to a can be written as 



2(7*2 



^ (re -ri) + Mi-Mc= ^ 



4o-J(rc -ri) 



G 



(10) 



where = 47r/9,;rj?/3, M, = Airpir^/S, r, = r\:(t = 0) 
is the radius of the (preexisting) core when the hardening 
phase starts at f = 0, and pi = p*{ri). The core radius 
then grows as 



rc{t) 



-^G(mi + m2) 



J{a) 

a{t) 



da. 



a 



(11) 



The binary separation quickly falls below Vc and subse- 
quent evolution is slowed down due to the declining stellar 
density, with a hardening time, 



th = \a/a\ = 



(12) 



that becomes increasingly long as the binary shrinks. The 
mass ejected increases approximately logarithmically with 
time, and the binary 'heats' background stars at radii 
Tc ^ a. If the hardening continues sufficiently far, gravita- 
tional radiation losses finally take over, and the two BHs 
coalesce in less than a Hubble time.^ In the scheme for 
core creation investigated here, at all redshift z < 5 the 
binary hardening timescale is always less than the time it 
takes the 'satellite' halo-|-BH system to sink to the center 
of the more massive progenitor due to dynamical friction 
against the dark matter background. 

The above relations assume that the stellar velocity dis- 
persion remains constant during the hardening of the bi- 
nary. This is a reasonable assumption as the two-body 
stellar relaxation timescale 



tr = 0.34- 



7Gyr(^ 
Vpc 



10 
inA 



CT150, 



is typically much longer than the hardening time. Here 
A is of order the total number of stars, m* is the typi- 
cal stellar mass, and cri5o is the stellar velocity dispersion 
in units of 150 kms~^. More importantly, we have ne- 
glected the depopulation of the loss cone, since it is the 
total stellar density that is allowed to decrease following 



equation (11), not the density of low angular momentum 
stars. The effect of loss-cone depletion (the depletion of 
low-angular momentum stars that get close enough to ex- 
tract energy from a hard binary) is one of the major un- 
certainties in computing the merger timescale, and makes 
it difficult to construct viable merger scenarios for BH bi- 
naries. The wandering of the binary center of mass from 
the galaxy center induced by continuous interactions with 
background stars (Quinlan & Hernquist 1997), the large 
supply of low-angular momentum stars in significantly flat- 
tened or triaxial galaxies (Yu 2002), the presence of a third 
BH (Vahonen et al. 1994; Blaes, Lee, & Socrates 2002; Pa- 
per I), the loss of orbital angular momentum to a gaseous 
disk (Gould fc Rix 2000; Armitage & Natarajan 2002), 
and the randomization of stellar orbits due to the infall 
of small satellites (Zhao, Haehnelt, & Rees 2001) will all 
mitigate to some extent the problems associated with loss- 
cone depletion (which may ultimately cause the binary to 
stall) and help the binary merge. 

4. RESULTS 

The ability of SMBH binaries in shaping the central 
structure of galaxies depends not only on the efhciency 
of three-body interactions at hardening the binary but 
also on how galaxy mergers affect the inner stellar den- 
sity profiles, i.e. on whether cores are preserved or steep 
cusps are regenerated during major mergers. To bracket 
the uncertainties and explore different scenarios we run 
three different sets of Monte Carlo realizations. In the 
first {'cusp regeneration') we assume, as in Paper I, that 
the stellar cusp oc is promptly regenerated after every 
major merger event, i.e., we replenish the mass displaced 
by the binary and reset = after every major merger. 
In the second ('core preservation^) the effect of the hier- 
archy of SMBH binary interactions is instead cumulative, 
i.e., ri is allowed to grow continuously during the cosmic 
evolution of the host. This second model is supported by 
N-body simulations involving mergers of spherical galaxy 
models with different density profiles, and showing that 
the remnant profile is quite close to the profile of the pro- 
genitors - in other words that the core appears to be pre- 
served during such mergers (e.g. Fulton & Barnes 2001). 
The third set is similar to the second, except that stellar 
ejection is switched off at separation a = 0.1 ah, as may be 
expected in the case the depletion of the loss cone were to 
'stair the binary around that separation (e.g. Milosavlje- 
vic & Merritt 2002). Thereafter the binary is assumed to 
shrink rapidly due to (say) gas processes and coalesce. As 
shown below, this case yields results that are intermediate 
between the first two. 

Let us focus on the first two sets of realizations. Fig- 
ure 1 shows how the mean core radius grows as a func- 
tion of redshift in the simulated merger history of a Mq = 
2 X 10^2 aj^fj a ]\^^ = iQia jjalo. The cores cre- 
ated by shrinking SMBH binaries tend to remain small 
(rc < 1 pc) until relatively recent times: on average, core 
radii approximately double between z = 1 and the present 
epoch. Core radii are also twice as large in the 'core preser- 
vation' case than in the 'cusp regeneration' one. (Note that 
in the latter averaging over all realizations smoothes out 



■^As in Paper I, the final mass rrtBH of the SMBH formed after coalescence assumes the entropy-area relation for BHs (maximally efficient 
radiative merging) : jj = + m'^. 
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the Ti = resetting after every major merger event.) 




0.01 ^ 



al. (1997). Our core sizes, however, are sensitive to the 
adopted density profile. We have assumed here that the 
continuous eroding action of shrinking binaries generates 
a flat density core. In the simulations of Milosavljevic & 
Merritt (2001), mergers of equal-mass stellar systems con- 
taining SMBHs and steep central density cusps produce 
nuclei with shallower cusps cx rather than flat cores. 



100 r 



Fig. 1. — Mean (averaged over 20 realizations) core radius of 
the main halo of the merger tree as a function of redshift. The 
expansion of the core is supported by the binding energy liberated 
by shrinking SMBH binaries. Solid line: Mq = 2 X lO^^ y;^ h^lo. 
Dashed line: Mq = 10^^ Mq halo. Left panel: cusp regeneration 
case. Right panel: core preservation case. 

Because of the lower central densities, binaries that form 
in the mergers of massive galaxies decay much less rapidly 
in the model that preserves cores than in the one where 
cusps are regenerated. Figure 2 shows the expected mean 
core radius at z = 0, as a function of halo mass. Small 
galaxies have tiny cores or no core at all, while massive 
galaxies have core radii that can exceed hundreds of par- 
sees. 

This is due to two effects: (1) our accretion recipe is 
fixed in order to reproduce the observed local tobh — f* 
relation. More massive halos have more massive central 
SMBHs and hence larger cores; and (2) the assumed 'bias' 
in the frequency of primordial seed BHs. Low mass galax- 
ies today are 'antibiased' as they arc assigned a smaller 
abundance of BH seeds with respect to the mean at early 
epochs, and some of them never host binaries. By con- 
trast, more massive galaxies form through the merging 
of an above-average number of primordial high-cr mini- 
halos, and so the damage done by shrinking BH bina- 
ries is enhanced. Note that, in the case of core preser- 
vation, the core radius is typically 3 times larger than 
the radius of the "sphere of influence" of the final SMBH, 
fBH — Gm^Yi/al. The predicted Tc of the most massive 
halos in the core preservation case is larger than the core 
one would simply compute from equation (11) assuming 
Ti = 0, an equal mass BH merger (toi = m2), and a final 
merged BH mass that satisfies the tobh ~ cr* relation of 
Ferrarese (2002), Tc oc tobh/o"^ oc niQ^ . 

How do our results compare with observations? The 
simple models for core creation described above yields 
core radii that scale almost linearly with galaxy mass, 

(X M°-^^°-^ in the range lO^^ < Mq < 4 x lO^^ 
A similarly scaling relation was observed by Faber et 

^More precisely, the mass deficit is defined as the difference in intej 
extrapolated inward from the break radius. 
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Fig. 2. — Mean core radius as a function of halo mass. Errorbars 
are l-cr rms. Lower dots: cusp regeneration case. Upper dots: core 
preservation case. Dashed line: core radius computed from eq. (11) 
assuming = 0, an equal mass BH binary (mi = m2), and a final 
merged BH mass that satisfies the rriBH ~ o"* relation. Note that 
this prescription overpredicts rc at low masses: this is due to small 
SMBHs deviating in our model from the extrapolated rriBH — c* 
relation. 

Observationally, even 'core' galaxies exhibit very shal- 
low power-law density profiles. A better test of our model 
predictions against galaxy data is provided by the 'mass 
deficit', i.e. the mass in stars that must be added to the 
observed cores to produce a stellar cusp.^ Milosavl- 
jevic et al. (2002) (see also Ravindranath, Ho, & Filip- 
penko 2002) have recently shown that the mass deficit in- 
ferred in a sample of early-type 'core' galaxies correlates 
well with the mass of their nuclear black holes, consistent 
with the prediction of coalescing SMBH binary models. In 
our model this quantity. 



Afdcf = 4(72r,/(3G), 



(14) 



is proportional to the total stellar mass ejected by shrink- 
ing binaries over cosmic history, and is directly related to 
the hardening efficiency. 

Figure 3 compares the mass deficit inferred from the 
data with the same quantity found at z = in our merger 
tree, 



Mdcf / Mq = (8.2 ± 3.8) (mBH/ Mq) 
(core preservation) and 



0.98±0.03 



(15) 



(11.5 ±4.2) (tobh/ M0)"-^3±o.O2 (16) 



Mdof / Mq 

;rated mass between the deprojected density profile and a r~^ 
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(cusp regeneration), where the fit was performed in the 
same men range of the observations. The cusp regener- 
ation case predicts Afdof — 2.7mBH for tobh = 10^ M©, 
and clearly underestimates the mass deficit observed in 
massive 'core' galaxies, (M^ct) — lOmeH- We find that 
typically about half of Afdcf is generated in the last binary 
merger. Note that the mass deficit defined in this way is 
about 25% larger than the stellar mass physically ejected 
in our simulations, since in our scheme ti* (hence the mass 
deficit) keeps increasing after binary coalescence as a con- 
sequence of the accretion of small dark matter satellites 
along the merger tree. 
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and switches off stellar ejection when a — 0.1 ah, i.e. stel- 
lar dynamical processes become inefficient at small binary 
separations. Even in this scenario stellar cusps may be 
efficiently destroyed over cosmic time by decaying SMBH 
binaries: the expected mass deficit is less than a factor 
of 2 smaller than in the standard core preservation case. 
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Fig. 4. — Same as Fig. 3, except that stellar ejection is switched 
off at separation a = 0.1 a^; thereafter the binary is assumed to 
shrink rapidly due to gas processes and coalesce. The realizations 
assume cores are preserved during galay mergers. 



Fig. 3. — Mass deficit produced at z = by shrinking SMBHs in 
our merger tree as a function of nuclear SMBH mass [filled dots). 
Left panel: cusp regeneration case. Right panel: core preservation 
case. Galaxies without a core (i.e. those that have never experienced 
a BH-BH merger or with their cusp recently regenerated) are shown 
as vertical arrows at an arbitrary mass deficit of 10*5 M0 . Empty 
squares: mass deficit inferred in a sample of galaxies by Milosavl- 
jevic et al. (2002). Empty circles: same for the 'core' galaxies of 
Faber et al. (1997) sample. Upper limits are shown for the 'cuspy' 
galaxies of Faber et al. (1997) sample, assuming a flat core within 
the upper limit on the core size. 

The predicted M^ef — wbh correlation is too tight com- 
pared to the observations: a distribution of inner den- 
sity profiles before and after galaxies merge and obser- 
vational errors may both contribute to the large scatter 
in the data points. We also draw attention to the fact 
that the mass deficits estimated for the cuspy galaxies of 
Faber et al. (1997) should be strictly considered as upper 
limits, and the observed M^cf — wbh relation may then 
be steeper than predicted (see also Ravindranath et al. 
2002). A few galaxies in the sample of Faber et al. have 
cores significantly smaller than expected from their cen- 
tral BH mass (e.g. NGC3115, with a break radius ~ 2 
pc and rTiBH = 10^ M©). Similarly the Milky Way galaxy 
has a BH with mBH = 3 x 10^ M© and a stellar core ra- 
dius Tc = 0.38 pc (Genzel et al. 2000). Within our simple 
scheme, objects with these properties are very rare in the 
'core preservation' case, and may require some form of 
cusp regeneration mechanism. 

The mass deficit measured in our third set of realizations 
is depicted in Figure 4. This assumes core preservation 



5. DISCUSSION 

We have assessed a hierarchical structure formation sce- 
nario in which galaxy inner cores are created from the 
binding energy liberated by shrinking SMBH binaries. The 
binary orbital decay heats the surrounding stars, erod- 
ing a preexisting stellar cusp oc r~^. A simple model in 
which the effect of the hierarchy of SMBH interactions is 
cumulative produces a correlation between 'mass deficit' 
and the mass of nuclear SMBHs with a normalization and 
slope that are comparable to the observed relation. This 
model is also able to reproduce the observed scaling re- 
lation between galaxy luminosity and core size. Relating 
the dark halo mass to galaxy luminosity using the mass-to- 
light ratios of van den Bosch, Mo, & Yang (2002), we find 
Vc cx Ly^ in the mass range 10^^ < Mq < 4 x 10^^ M©, in 
agreement with the scaling observed by Faber et al. (1997). 
Despite these encouraging results, models of core forma- 
tion by binary SMBHs remain uncertain. The ability of 
SMBH binaries in modifying the inner density profiles of 
galaxies depends on whether cores are preserved or steep 
cusps are regenerated during galaxy mergers, and on the 
poorly known efficiency of three-body interactions at hard- 
ening the binary. 

There are of course further complications: (1) If the age 
of the stellar core-FBH system is much larger than the local 
relaxation time, then an equilibrium distribution of bound 
stars will be set up within the BH sphere of influence. We 
find that such a steady-state density cusp, p» oc r^^/^ 
(Peebles 1972; Bahcall & Wolf 1976; Frank & Rees 1976; 
Lightman & Shapiro 1977) will have no time to develop 
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in our hierarchical scheme for the assembly of SMBHs, as 
the relaxation timescales are typically rather long. If the 
age of the system is instead much less than the relaxation 
time, and the hole grows by accreting gas on a time scale 
long compared with the orbital period of the surrounding 
stars, then again a power-law density cusp will be cre- 
ated, this time with slope r~^/^ extending to Tbh (Young 
1980). As the mass of the hole increases, the core will start 
to loose its identity until, when tbh ^ fc the cusp will 
join smoothly onto the isothermal profile. Adiabatic 
growth of stellar cusps is typically not important within 
our scheme, since we find that the core mass Mc exceeds 
mBH, i-e. that tbh < fc in most systems; (2) Dark matter 
particles will be ejected by shrinking SMBH binaries in 
the same way as the stars, i.e. through the gravitational 
slingshot. Their contribution to binary hardening will be 
weighted by the fraction /dm = /9dm/(/Cdm+P*), typically 
<C 1 in galaxy cores for our assumed density distributions. 
The NFW profile is shallower than the SIS profile in the 
inner regions, so ejection of a given mass will create a 
larger core in an NFW profile. The mass ejected in dark 



matter is, however, significantly less than the mass ejected 
in stars. Eventually, the destruction of dark matter cusps 
by binary SMBHs may lead to cores of the same extent 
as the stellar ones. Furthermore, the erosion of the stellar 
cusp may lower the central dark matter density through a 
process ('adiabatic expansion') that is the opposite of the 
adiabatic compression of dark halos in response to bary- 
onic cooling and infall (Blumenthal et al. 1986). While tlu^ 
observed rotation curves of dwarf and lower surface bright - 
ness galaxies suggest that the inner regions have constant 
density cores rather than the density cusps predicted by 
CDM'(McGaugh & De Blok 1998; Swatcrs et al. 2003), 
shrinking BH binaries will not alleviate significantly this 
'CDM cusp problem', since the predicted mass deficit is 
much lower than that required to fit the rotation curves. 

We have benefitted from discussions with G. Blumen- 
thal, S. Faber, M. Haehnelt, M. Milosavljevic, and M. 
Recs. Support for this work was provided by NASA 
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0205738 (P.M.). 
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